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Osmosis plays a central role in the function of living and soft matter systems. While the thermo- 
dynamics of osmosis is well understood, the underlying microscopic dynamical mechanisms remain 

o, 

the subject of discussion. Unravelling these mechanisms is a crucial prerequisite for eventually 
understanding osmosis in non-equilibrium systems. Here, we investigate the microscopic basis of 
§ '. osmosis, in a system at equilibrium, using molecular dynamics simulations of a minimal model 

in which repulsive solute and solvent particles differ only in their interactions with an external 

S : 

potential. For this system, we can derive a simple virial-like relation for the osmotic pressure. Our 
simulations support an intuitive picture in which the solvent concentration gradient, at osmotic 
equilibrium, arises from the balance between an outward force, caused by the increased total den- 
sity in the solution, and an inward diffusive flux caused by the decreased solvent density in the 
solution. While more complex effects may occur in other osmotic systems, they are not required 



for a description of the basic physics of osmosis in this minimal model. 
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I. INTRODUCTION 



Osmosis is the physical phenomenon in which a concentration difference of impermeant 
solute molecules across a semi-permeable membrane produces a difference in solvent density, 
and in pressure, across the membrane. This has enormously important consequences - 
osmosis forms the basis of the transport of ions across cell membranes l| , the regulation of 



blood salt levels by the kidneys |2j, uptake of water by plants |3|, |4|, technologies for kidney 
dia, ysis B and clean power paction gg and many other p_. Osmosis a,s„ lies at 
the heart of many important phenomena in chemical physics including the Donnan effect [7[, 
the depletion interaction |8[ and, recently, mechanisms for generating self-propelled colloidal 
particles |9l-|l5|. 

The cornerstone of our understanding of osmosis is the van 't Hoff relation, which states 
that a solute concentration difference, Ac s , leads to an equilibrium osmotic pressure differ- 



ence AP = ksTAcg 



16 



171 ] . Remarkably, this relation predicts that the osmotic pressure 



difference is the same as the pressure of an ideal gas at concentration Ac s , regardless of the 
molecular nature of the solute and solvent molecules. Gibbs showed that the van 't Hoff 
relation can be derived by setting the chemical potential of the solvent equal across the 
membrane 
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2l| ; the thermodynamic origin of the osmotic pressure difference is then the 



entropy of mixing between solute and solvent. The van 't Hoff relation only holds for low 
solute concentrations; at higher concentrations, deviations from this "ideal" behaviour can 

n 

be used as a diagnostic of solute-solute interactions, such as ion pairing [22]. 



Thus, from a thermodynamic point of view, osmosis is wel 



understood. However, osmosis 



23|, for which a thermodynamic 



is also important in non-equilibrium systems 
description is not valid. A clear picture of the underlying molecular dynamical mechanisms 
is an essential prerequisite for understanding osmotic phenomena in systems that are out 
of equilibrium. Yet such a picture is largely lacking, even though the dynamical basis of 
osmosis has been the subject of over 100 years of discussion, since the seminal work of van 
't Hoff. Relevant factors may include bombardment of the membrane by solute molecules, 
diffusion of solvent across the membrane driven by its density difference, specific types of 



interaction between solute and solvent molecules and pu 



membrane in the wake of solute-membrane collisions 



l ing of solvent molecules across the 



2im\- 



In an important contribution, Brady 34j showed that for a system of colloidal solutes 
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suspended in a coarse-grained solvent, the osmotic pressure can be derived from a purely 
mechanical, hydrodynamic Langevin description of colloidal motion; here, hydrodynamic 
interactions between colloids, mediated by solvent, play a central role. For systems with 



explicit solvent, Guell and Brenner 



35| showed that the van 't Hoff relation can be de- 



rived from first principles by taking into account all the forces involved; here, the role of 
membrane-solute interactions is emphasized. Importantly, while a thermodynamical de- 
scription of osmosis cannot depend on the details of the system's dynamical rules (e.g. the 
interparticle interactions, membrane geometry, etc), the nature of the underlying dynamical 
mechanisms may be sensitive to these details, and could be different for different osmotic 
systems. Indeed, detailed models of osmotic water transport through membrane pores in 
ionic solutions have highlighted the importance of specific solvent-membrane and solvent-ion 
interactions in these systems 36M38| , and the unusual characteristics of osmotic water flow 



through carbon nanotube membranes, which can be described by ID stochastic hopping 
models [39(] . Other studies have shown that partial solute penetration into membrane pores 
can strongly influence solvent flow 40, |4J|, and that, for ionic solutes, non-trivial coupling 
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44|. 



can arise between the transmembrane electric field and osmotic flow 

In this paper, we seek to provide a clear picture of the dynamical mechanisms underlying 
osmosis, in a system which is as simple as possible. Murad and Powles 
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47] . and Itano 



et al [48|, have demonstrated that osmosis can be achieved in systems in which solute and 
solvent particles have identical Lennard- Jones interactions. Here, using an even simpler 
system with purely repulsive particles, and with no solvent-membrane interactions, we carry 
out a detailed study of the dynamics of solvent and solute particles in the osmotic steady 
state. For this system, we are able to derive a simple expression for the osmotic pressure 
in terms of solute-solute and solute-solvent interactions. Our results point to a simple and 
intuitive model in which the solvent density and pressure gradients in osmotic equilibrium 
are maintained by a balance between an outward solvent flux generated by the osmotic 
pressure gradient and an inward diffusive flux. Although more complex mechanisms may be 
at play in real osmotic systems, they are not required to describe the fundamental physics 
of osmosis in our system. We hope that this work will provide a basis on which to build 
future descriptions of osmotic phenomena in both equilibrium and non-equilibrium systems. 
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II. A MINIMAL MODEL FOR OSMOSIS 



Our model system is designed to reproduce the essential physics of osmosis, while elimi- 
nating as far as possible all complications due to the specific chemical nature of the solutes, 
solvent and the semi-permeable membrane. In our system solute and solvent particles inter- 
act with each other indistinguishably, and there are no solvent-membrane interactions. Our 
simulation setup is illustrated in Figured] The simulation box is partitioned into "solution" 
and "solvent" compartments; the solute particles are confined within the central solution 
compartment by an external potential of the form U = ksT{a / dx)~ 9 (where dx is the per- 
pendicular distance between a solute particle and the boundary of the compartment, and a 
is the particle diameter) 59[. The solvent particles do not experience the confining poten- 
tial and are free to move throughout the simulation box. The confining potential therefore 
acts like a semi-permeable membrane. Both solvent and solute particles are of unit mass 
m and interact via identical, repulsive, Weeks, Chandler, Andersen (WCA) interactions 



(U(r) = 4e 
e = k B T = 



(jr) 12 — (^) 6 + \ if r < 2 1 / 6 cr and zero otherwise), with parameters o = 1 and 
49[ (i.e. our units of energy and distance are k B T and a respectively^ 



system is simulated using Molecular Dynamics, with the velocity Verlet algorithm 50 



'he 



and timestep 0.001 (in reduced units 50| 60j), combined with a Nose- Hoover thermostat 



50 



521 ] . In total, our system contains 5000 particles (solute plus solvent) in a periodic, cubic 
simulation box of size L = 18.42cr, so that the total particle density ptotai = 0.8cr~ 3 (the 
packing fraction np tota \/6 = 0.42). 

In this paper, we focus on the osmotic steady state; we therefore allow the system to 
equilibrate thoroughly before data is collected [61(- We measure the local pressure in the 



solution and solvent compartments using the Method of Planes [53|, [54j. As discussed in 
Appendix |A[ this method constitutes a direct measurement of the kinetic and interaction 
components of the momentum flux across a local plane. To compute the concentration of 
solute particles, c s , we need to define the volume of the solution compartment. This is 
non-trivial, because the confining potential is smooth. We define the solution volume by 
matching the pressure-density relation for a "gas" of solute particles, confined in the solution 
compartment, in the absence of solvent, with that of a system of WCA particles simulated 
in a periodic box (see Appendix [A]) . 
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FIG. 1: Illustration of our model system, (a): Solute particles are confined within a cubic 

solution compartment located in the centre of the simulation box. (b): Simulation 
snapshot; solute and solvent particles are coloured green and blue respectively. For clarity 
only particles located in the solution compartment are shown. 

III. OSMOSIS IN THE MODEL SYSTEM 

A. Density imbalance 




la x /a 



(a) (b) 

FIG. 2: Local density profiles p(x) (in units of er~ 3 ) measured across the middle of the 
simulation box, for solute concentration c s = 0.254cr~ 3 . Panel (a) shows the total particle 
density, ptotai^); panel (b) shows the local solvent density p v (x). 

We first verify that osmosis indeed occurs in our model system by measuring the steady- 
state local density profiles in the solution and solvent compartments. Figure |2] shows density 
profiles, taken through the middle of our simulation box. As expected, the total particle 
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FIG. 3: Spatially averaged solvent density p v as a function of solute concentration c s (both 
in units of er~ 3 ), in the solution and solvent compartments (circles and squares 
respectively). The dashed lines show theoretical predictions based on the 
Carnahan-Starling equation of state (see Appendix [Bj . 

density, solute plus solvent (Figured), is higher in the solution than in the solvent com- 
partment. In contrast, the solvent density is lower in the solution than in the solvent 
compartment (Figure Eb). This density imbalance increases linearly with the concentration 
of solute, Cg, as shown in Figured These results are in good agreement with thermodynamic 
predictions, obtained by approximating the WCA particles as hard spheres and setting the 
solvent chemical potential, obtained from the Carnahan-Starling equation of state, equal in 
the two compartments (dashed lines in Figure O for details of the calculations see Appendix 
IB]) . Figure |3] also demonstrates that over the parameter range considered here, the particle 
density in the solvent compartment remains virtually unaffected by changes in the solute 
concentration, confirming that it is large enough to be regarded as a reservoir. 

B. Osmotic pressure 

The osmotic pressure, AP (i.e. the pressure difference between the solution and solvent 
compartments), is shown in Figure H] as a function of the solute concentration, c s . At low 
solute concentration, our simulation results (circles in Figure H|) are in good agreement 
with the van 't Hoff relation (dotted line). At high solute concentration (c s > 0.1o~~ 3 , 



0.8 



> 
Q_ 



0.6 
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corresponding to solute packing fraction greater than 0.05), the osmotic pressure exceeds 
that predicted by the van 't Hoff relation, as expected in a system with repulsive solute- 
solute and solute-solvent interactions. Interestingly, however, the osmotic pressure in our 
system is significantly lower than the pressure that would be obtained for a gas of WCA 
particles at density c s (dashed line in Figure H]). Thus a naive picture in which one treats 
the solution simply as a "solute gas", ignoring the solvent, does not correctly account for 
the osmotic pressure. 




FIG. 4: Osmotic pressure AP (in units of ksTa~ 3 ) as a function of solute concentration c s 
(in units of <r~ 3 ). In the main plot, the circles show simulation results, computed using the 
method of planes (see Appendix [A}, and the squares show AP computed using Eq. (CQ). 
The dotted line shows the van 't Hoff relation, AP = fc^Tc s , while the dashed line shows 

the pressure of an equivalent system of hard spheres at density c s , computed using the 
Carnahan-Starling equation of state. In the inset, the circles are the same as in the main 

plot, while the triangles show the prediction of our simple hopping model (Eq.(j2])) 
AP = —kBTp° nt log [/?v n /P? ut ] (where and p° ut are the average solvent densities in the 
solution and solvent compartments, respectively). 



We show in Appendix O that, starting from the Clausius virial relation for the solute 
particles, an expression for the osmotic pressure can be derived in terms of the forces of 
interaction between the solute and solvent particles: 



AP » k B Tc s + + E E ) 

\ is js>is is jv / 



(1) 



Here, V is the volume of the solution compartment (defined as described above), fi is the 
position of particle i and /y is the force exerted on particle i by particle j. The first term 
in the brackets sums over all pairs of solute particles (avoiding double counting) while the 
second term sums over all solute-solvent pairs (note that here particle i denotes the solute 
and j the solvent; the solvent may be inside or outside the solution compartment). Since the 
starting point of this derivation, the Clausius virial relation, follows from Newton's equations 
of motion in a system at steady state (see Appendix [C]), Eq. ([1]) amounts to a mechanical 
description of the osmotic pressure in our system. Figure 0] (squares) shows that Eq. ([I]) is 
indeed in good agreement with our previous computation of the osmotic pressure (circles), 
over the full range of solute concentrations tested. 




FIG. 5: Contributions to the deviation of the osmotic pressure from the van 't Hoff 
relation. The circles show the solute-solute contribution (first term in the brackets in 
Eq. (pQ)); the squares show the solute-solvent contribution (second term in the brackets in 
Eq. ([!])). The inset splits the solute-solvent contribution into the parts due to interactions 
with solvent particles inside and outside the solution compartment (upward and 
downward-pointing triangles, respectively). 



Eq. ([I]) provides insight into the origins of the deviation of the osmotic pressure from 
the van 't Hoff relation as the solute concentration increases. This deviation arises from 
both solute-solute and solute-solvent interactions (the two terms in brackets in Eq. ([[]), 
respectively). These contributions are plotted in Figure |5j the solute-solute interactions 
(circles in Figure|5]) make a positive contribution to the osmotic pressure, as one would expect 
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from a "solute gas" picture, with repulsive solute-solute interactions. However, solute-solvent 
interactions (squares in Figure |5]) make a net negative contribution, decreasing the osmotic 
pressure. To understand the origin of this negative contribution, we split the solute-solvent 
term into contributions due to solvent particles inside and outside the solution compartment. 
The inset to Figure shows that solute-solvent interactions inside the solution compartment 
increase the osmotic pressure (upward triangles in the inset), while interactions between 
solute particles and solvent particles which are outside the solution compartment (i.e. cross- 
membrane interactions; downward triangles in the inset) decrease the osmotic pressure. 
Thus, in our system, solute particles are pushed outwards by interactions with other solutes 
and with solvent particles in the solution, but at the same time they are pushed inwards by 
interactions with solvent particles located outside the solution compartment. For colloidal 
dispersions, in which the solutes are colloidal particles, it is usual to treat the system as a 
single- component fluid of colloids which interact via an effective potential that takes into 
account the effects of the solvent. Because the colloids are orders of magnitude larger 
than the solvent molecules, the effective interactions are well-represented by a pairwise 
intercolloidal potential V(r). By analogy with atomic systems the osmotic pressure, IT, is 
then written as II = nksT — ^-n 2 J °° r 3 g(r) (dV/dr) dr, where n is the number density of 
colloidal particles and g(r) is the radial pair distribution function BQQ. This expression 
can a. 8 o be derived from a purely dynauncal description of colloidal motion, taking account 
of Brownian motion and hydrodynamic interactions [341 ]. Eq. (CLJ) provides an analogous 
expression for the osmotic pressure, for a system in which the solvent degrees of freedom are 
treated explicitly and on the same footing as those of the solute. 



IV. WHAT MAINTAINS THE SOLVENT DENSITY GRADIENT? 



A. Balance between outward and inward fluxes 



Our simulations allow us to investigate in detail the microscopic forces which produce the 
imbalance in the density of solvent particles between the solution and solvent compartments. 
Figure [6] shows the net force per particle, acting on the solvent particles, as a function of 
position x, in a slab through the middle of the simulation box. Solvent particles close to the 
boundaries of the solution compartment are on average pushed out of the solution, towards 
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FIG. 6: Net force per particle, normal to the membrane (in units of /c^Ter -1 ) acting on the 
solvent particles, as a function of position x, in a slab taken through the middle of the 
simulation box, for solute concentration c s = 0.254a -3 . Solvent particles close to the 
boundaries of the solution compartment tend to be pushed out of the solution. 



the solvent compartment. This net force arises from the fact that the total density inside the 
solution compartment is higher than that in the solvent compartment (as shown in Figure 
[2k) - thus, we expect solvent particles at the boundary to experience more collisions from 
the solution side than from the solvent. Why then is there no net flow of solvent out of the 
solution compartment in response to the net outward force? 

The answer lies in the solvent density profiles. As shown in Figure [2b, the density of 
solvent particles is higher in the solvent compartment than in the solution. This creates a 
diffusive flux of solvent particles into the solution, which counteracts the outward flux caused 
by the net outward force. From a simplistic point of view, let us consider the "hopping" 
of solvent particles across the boundary. An individual particle on the solution side has a 
high probability to leave the solution, due to the outward pushing force, while an individual 
particle on the solvent side has only a low probability to enter the solution. However, there 
are more solvent particles on the solvent side and consequently more attempts to hop into 
the solution than out of it. Thus the net fluxes of solvent particles inward and outward are 



equal and the solvent density imbalance is maintained 62]. 

Using this picture, we can relate the solvent density imbalance to the osmotic pressure. 
Let us assume the solvent densities, and p° ut in the solution and solvent compartments 
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respectively, are uniform, with a very sharp density change at the boundary. We now define 
slabs of width h (of the order of the molecular size) next to the boundary on the solvent and 
solution sides, and consider the hopping of solvent particles between these slabs. Solvent 
particles that hop from the solvent side into the solution experience an energy penalty hf out , 
where the outward force / out arises from the fact that the total density is higher in the 
solution, and can be related to the osmotic pressure AP by AP = / ou t/a, where a is the 
slab area per particle. Setting the fluxes of solvent particle "hops" equal in the inward and 
outward directions, we find that = p° ut exp -/ ou tV(^B^) = P° ut ex P [ — ahAP / {ksT)]. 
Noting now that ah is the volume per particle, which is approximately the same on both 
sides of the membrane - i.e. that ah ps 1/Ptotai ~ l/p° 



?° ut , leads to the result 



AP 



P 



.out 



(2) 



-k B TpT \og 

Eq. (j2J) seems to make sense from a thermodynamic point of view, since it amounts to 
a Boltzmann relation p™/p° ut = exp[-AU/k B T], with AU = AP/p° ut . Figure H (inset; 
triangles) shows that Eq. (TS]) holds in our simulations, even for solute concentrations well 
beyond the limit of validity of the van 't Hoff relation. For low solute concentrations, Eq.(J2]) 
reduces to the van 't Hoff relation. Noting that = p^otai — c s ~ P° ut ~ c s leads to: 



log 



^out 
rv J 



AP 



(3) 



k B Tp^ 

The van 't Hoff relation is then recovered when we expand the logarithm to first order in 
c s /p° ut (which is small for low solute concentration): 



AP^T/CMog 



P 



.out 



/CrTCs. 



(4) 



B. Solvent-solute correlations 



The simple picture presented above provides a kinetic description which is consistent 
with the known thermodynamic results, and which is also consistent with the force profiles 
measured in our simulations. To show that this actually describes what is happening in our 
simulations, we need to rule out more complex dynamical mechanisms. A key feature of our 
simple picture is that it does not predict any correlations between the configuration of the 
solute particles and the "hops" of solvent particles across the membrane - unlike some other 
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proposed mechanisms in which solvent particles are "pulled" across the membrane in the 
wake of nearby solute particles. To look for such effects in our simulations, we therefore test 
whether correlations between solute and solvent motion influence the dynamics of solvent 
particles close to the membrane. 

We first ask whether the microscopic flux of solvent particles into the solution compart- 
ment is influenced by the presence of nearby solute particles. At time i, we focus on solvent 



63| outside the solution com- 



particles which are situated less than 0.001 particle diameters 
partment. Within the next timestep, some of these solvent particles cross into the solution 
compartment, while others do not. Figure [Tti shows the distribution of distances to the 
nearest solute particle, for those solvent particles which do cross in the next timestep (black 
line) and for those which do not (red line). These distributions are indistinguishable, sug- 
gesting that proximity of a solute particle does not make any difference to the chance that 
a solvent particle, located close to the boundary, will cross into the solution compartment. 

0.8F i 

0.4- 





FIG. 7: Probability distributions for the distance d to the nearest solute particle (a) and 
the perpendicular velocity v of the nearest solute particle (b), for solvent particles located 
less than 0.001 particle diameters outside the solution compartment, which do (black lines) 
or do not (red lines) enter the solution compartment in the next simulation timestep. In 
these simulations, c s = 0.05<7~ 3 . The boundary of the solution compartment is defined as 

described in Appendix |A] 



One might also suppose that the velocity of nearby solute molecules might influence the 
microscopic flux of solvent particles - e.g. a solvent particle could cross into the solution in 
the wake of a nearby solute particle as it rebounds from the membrane. Figure [7b shows the 
probability distribution for the velocity of the nearest solute particle (perpendicular to the 
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solution boundary), for solvent particles located close to the solution boundary, which do and 
do not subsequently cross into the solution. Once again, the two probability distributions 
are indistinguishable. 



0.8 




3 



FIG. 8: (a): 9 is defined as the angle between the instantaneous velocity of a solvent 
particle i as it enters the solution compartment and the vector joining it to particle j 
within the solution compartment, (b): Probability distribution P(9) when particle j is 
defined to be the nearest solute particle (black line) and when particle j is the nearest 
particle of either solute or solvent (red line). In these simulations c s = 0.05cr~ 3 . The 
boundary of the solution compartment is defined as described in Appendix [A] 

As a final test, we measured whether solvent crossings occur preferentially in the direction 
of the nearest solute particle. If solute particles tend to "suck" solvent particles across the 
membrane, solvent crossings should be preferentially oriented towards the nearest solute 
particle. We therefore measured the angle 9 between the instantaneous velocity Vi, of a 
solvent particle as it crossed into the solution compartment, and the vector fji = fj — f\, 
where and fj are the position of the solvent particle and the nearest solute, respectively 
(see Figure Et). Figure [8b (black line) shows the probability distribution P{9). We then 
repeated the calculation, but this time allowing particle j to be either solute or solvent 
(within the solution compartment). The resulting distribution is given by the red line in 
Figure [8b. Again, the two distributions are indistinguishable, showing that solvent particles 
do not cross the membrane preferentially in the direction of nearby solute particles. We 
therefore conclude that correlations between solute and solvent dynamics do not play a 
significant role in determining the microscopic solvent flux in our simulations. This suggests 
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that the simple dynamical picture described above is indeed sufficient to describe osmosis 
in our system. 



V. DISCUSSION 

In this paper, we have used a minimal model of an osmotic system to investigate the 
microscopic dynamical mechanisms underlying the maintenance of equilibrium osmotic pres- 
sure and density gradients. In our system, solute and solvent particles interact identically, 
and there are no solvent-membrane interactions. For this system, we are able to derive 
an expression (Eq.flT])) for the osmotic pressure, in terms of the forces of interaction be- 
tween solute-solute and solute-solvent pairs of particles; this derivation does not assume any 
separation of scales between the solvent and solute particles. We find that, in our simula- 
tions, solute-solute and solute-solvent interactions are approximately equally important in 
determining the osmotic pressure - in contrast with the naive view that an osmotic system 



behaves like a "solute gas" , independent of the solvent [64 ] . 

Analysis of the density and force profiles in our simulations leads us to propose a very 
simple picture, in which the density imbalance of solvent particles across the membrane 
is maintained by a balance between an outward force-driven flux (due to the higher total 
density in the solution) and an inward diffusive flux (due to the lower solvent density in the 
solution). We show that a simple calculation, based on the "hopping" of solvent molecules 
across the membrane, leads to a relation between the solvent density gradient and the 
osmotic pressure (Eq.©), which is in good agreement with our simulation results even for 
high solute concentrations. For low solute concentrations, this expression reduces to the van 
't Hoff relation. We do not detect any evidence of correlations between the solute particles 
and the solvent flux, leading us to believe that this simple "hopping" picture is sufficient to 
describe the dynamical basis of osmosis in our simulations. 

In previous work, various molecular mechanisms have been implicated in osmosis, in- 
cluding solute-membrane collisions, diffusion of solvent across the membrane, solute-solvent 
interactions and correlations between solvent and solute dynamics. In simulations of realistic 



osmotic systems, the dynamics of solvent molecules inside pores 39|, |56|, interactions with 



the pore wall 



361 ] . interactions between solvent and solute molecules 



36 



and correla- 



tions between solvent flux and solute dynamics [40j have all been shown to play important 
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roles. It is possible that different dynamical mechanisms are at play in different osmotic 
systems, all of these being consistent with the known, and universal, thermodynamic rela- 
tions. Thus, an important topic for further work will be to understand how the conclusions 
of our work change for more realistic scenarios: e.g. when the membrane presents a barrier 
also to the solvent molecules (via explicit pores or via a smooth potential barrier), when the 
solute and solvent do not interact identically, etc. Our aim in this work, however, was to 
focus on a minimal model system, for which we could characterize in detail the underlying 
molecular mechanisms. For this system, we find that complex mechanisms are not required 
to understand the basis physics of osmosis. 

We hope that this study will provide a basis for developing our understanding of osmosis 
in systems out of equilibrium, for which thermodynamic concepts are not available. Recently, 
much attention has focused on active particles such as motile bacteria, swimming "Janus" 



colloids or colloidal "chuckers" 



151 ] . These can act as solutes: an interesting question 



concerns how the osmotic pressure of a suspension of active colloids differs from that of 
the corresponding passive suspension. Moreover, active particles may themselves generate 
osmotic gradients, with interesting potential for motility and self-organisation. The approach 
described here can easily be extended to investigate the dynamics of osmotic systems in which 
either the solute or the solvent particles are active. 

Finally, we note that the minimal setup presented in this paper also provides a very 
general way to measure osmotic pressure in systems involving solvent and solutes, simply by 
partitioning a system into two compartments using an external potential permeable only to 
the solvent. This may provide an alternative route, for example, to the detection of phase 
transitions, in cases where the traditional thermodynamic approaches are not feasible. 
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Appendix A: Pressure measurement and solution boundary 
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FIG. 9: (a): Method of planes with a plane of finite extent in the y and z directions. Solid 
unidirectional arrows denote contributions to the kinetic part of Eq. flAll) (</>); solid 
bidirectional arrows denote contributions to the interaction part of Eq. llAlj) . All 
interactions in which the line connecting the two particles crosses the plane are included 

(e.g. the bottom pair of particles); however the interaction between the top pair of 
particles (dotted arrow) is excluded, (b): Pressure-density relation of a "gas" of solute 
particles, confined in the solution compartment in the absence of solvent (closed circles), 
compared with that of a periodic box of equivalent particles (open circles). The pressure is 
in units of ksTa~ 3 ; the density is in units of cr~ 3 . For the closed circles, the density of the 
solute gas is defined using an effective boundary of the solution compartment positioned 
0.91a inside the locus of divergence of the the confining potential. 

We measure the local pressure in the solution and solvent compartments using the method 
of planes 0, 3]. In this method, one defines a plane within the region of interest, and 
monitors the time-averaged flux cf) of particle momentum normal to the plane (due to particles 
crossing the plane), and the interparticle forces normal to the plane, for pairs of particles 
on opposite sides of the plane. For a plane at position x = x', which extends right across 
the simulation box in the y and z directions, the pressure P is given by 




(Al) 



i,j>i 
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where A is the area of the plane, the sum is over all pairs of particles (without double 
counting), Xi and Xj are the x-components of the particle position (x being the normal to 
the plane) and f X ij is the normal component of the force between particles % and j. The term 
in the square brackets ensures that only pairs of particles which are on opposite sides of the 
plane contribute to P. In our simulations, we modify this method in order to use planes 
with finite extent in the y and z directions. In this case, the sum in Equation IA1I should 
include only those pairs of particles for which the line connecting the particle positions 
crosses the plane (see Figure [9^)- We have verified that this method produces results in 
good agreement with global pressure measurements (for homogeneous systems) and with 



other methods for measuring the local pressure [57J, even when the plane is very small, of 
the order of a few particle diameters. Interestingly, in our system, since there are no solvent- 
membrane interactions, the osmotic pressure is also given by the magnitude of the forces of 
confinement on the solute particles, per unit area of the confining boundary. This provides 
an alternative way to measure the osmotic pressure, with results in good agreement with 
those of the method of planes. 

To define the concentration of solute particles in our simulations, we need to define an 
effective boundary for the solution compartment. To do this, we match the pressure- density 
relation of a "gas" of solute particles, confined in the solution compartment in the absence 
of solvent (shown in Figure |9b, closed circles), with that of a periodic box of equivalent 
particles (Figure |9b, open circles). This results in an effective boundary of the solution 
compartment which lies 0.91<r inside the locus of divergence of the the confining potential. 

Appendix B: Hard sphere prediction for density imbalance 

In Figure El we compare our simulation results for the solvent density imbalance between 
solution and solvent compartments, to thermodynamic predictions for a system of hard 
spheres, obtained u sing the Carnahan-Starling expression for the free energy of a system 



of hard spheres 55|, |58). The solvent compartment, of volume V out , contains iV° ut solvent 
particles, at packing fraction r/° ut (this is related to the density p° ut by r/° ut = np° nt a 3 /6). 
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The free energy is then 



k R T 



N out 
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where A is the de Broglie thermal wavelength. The chemical potential is given by /i° ut 
dF ont /dN° nt = (7icr 3 /6V out )dF ont /dri° ut : 

"8 - 9r/° ut + 3r/° ut 



.out 

k B T 



log ( P e ut A 3 ) + V°, 



I-77: 



out\ 3 
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The solution compartment, of volume V m , contains a mixture of solute particles, at packing 
fraction r/ s , and solvent particles, at packing fraction 77™. Its free energy is 
pin 



k R T 



N s [log (p s A 3 ) - 1] + AC [log (/^A 3 ) - 1] 
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The chemical potential /< n = dF' m /dN™ = (na 3 /6V in )dF in /dr]^: 
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Setting /i™ = /x° , we arrive at the following relation: 



log 



out 



8 - 9(r/ s + r/i n ) + 3(r] s + rfc) 



out 



(l-(r/ s + C)) 3 
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Conservation of solvent particle number implies that 



(B5) 



out out 



V m rf* + V out ri: 
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(B6) 



where N v is the total number of solvent particles. Eqs. ( 1B5j) and (1B6[) are solved numerically 
to obtain 77™ and r/° ut as functions of the solute packing fraction 77 s , as plotted in Figure [3) 



Appendix C: Derivation of Eq.flT]) for the osmotic pressure 

Eq. ([1]) follows from the dynamical equations of motion for the solute particles, under 
the conditions that the system is in steady state and the solute particles are confined in 
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the solution compartment. The first of these conditions implies that d(^2 i fl ■ p^j jdt = 
(where the sum is over all solute particles and and Pi are the position and momentum 
of particle i respectively). Expanding the time derivative and using Newton's equations of 
motion leads to: 

Ef) + (X>^° (CD 

where Fj is the total force acting on solute particle i, and rru is its mass. Eq. (IClj) constitutes 



a Clausius virial relation for the solute particles [50j, |55|, |57|. Splitting the force into the 
contributions due to interactions with the confining potential, solute-solute interactions and 
solute-solvent interactions, and using the equipartition theorem, we arrive at: 

3N s k B T + ^fn + Yl f *- >w = °> ( C2 ) 

is js>is is jv is 

where is the force exerted on particle % by particle j and /i )Ccm f is the confining force on 
particle i. The first sum is over solute-solute pairs (without double counting), the second sum 
is over solute-solvent pairs, and the final sum is over interactions between solute particles 
and the "membrane". Focusing on the final term, we note that, to a good approximation, 
solute particles only feel the confining force when they are very close to the boundary of 
the solution compartment. The final term can then be approximated by an integral over 
a narrow shell at the boundary of the solution compartment. Taking into account also the 
fact that the osmotic pressure is given by the average normal force per unit area on the 
membrane, which in our system is equivalent to the magnitude of the confining force on the 



solute particles, per unit area, leads to 



X^/Unf ^3^ in AP (C3) 

is 

where AP is the osmotic pressure and V m is the volume of the solution compartment. 
Substituting Eq. (1C3jl into Eq. (1C2l) leads directly to Eq.flTJ in the main text. 
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